Astronomy <fe Astrophysics manuscript no. AA8495-accepted 


© ESO 2008 


February 2, 2008 





Leptonic secondary emission in a hadronic microquasar model 

M. Orellana^'^'*, P. Bordas^, V. Bosch- Ramon^ , G. E. Romero^'^'**, and J. M. Paredes^ 

^ Facultad de Ciencias Astronomicas y Geofi'sicas, Universidad Nacional de La Plata, Paseo del Bosque, 1900 La Plata, Argentina 

Institute Argentine de Radioastronomi'a, C.C.5, (1894) Villa Elisa, Buenos Aires, Argentina 
^ Departament d'Astronomia i Meteorologia, Universitat de Barcelona, Marti i Franques 1, 08028, Barcelona, Spain 
* Max Planck Institut fur Kernphysik, Saupfercheckweg 1, Heidelberg 69117, Germany 

Accepted 04/09/2007 

ABSTRACT 

Context. It has been proposed that the origin of the very high-energy photons emitted from high-mass X-ray binaries with 
jet-like features, so-called microquasars (MQs), is related to hadronic interactions between relativistic protons in the jet and 
cold protons of the stellar wind. Leptonic secondary emission should be calculated in a complete hadronic model that include 
the effects of pairs from charged pion decays inside the jets and the emission from pairs generated by gamma-ray absorption in 
the photosphere of the system. 

Aims. We aim at predicting the broadband spectrum from a general hadronic microquasar model, taking into account the 
emission from secondaries created by charged pion decay inside the jet. 

Methods. The particle energy distribution for secondary leptons injected along the jets is consistently derived taking the energy 
losses into account. We also compute the spectral energy distribution resulting from these leptons is calculated after assuming 
different values of the magnetic field inside the jets. The spectrum of the gamma-rays produced by neutral pion-decay and 
processed by electromagnetic cascades under the stellar photon field. 

Results. We show that the secondary emission can dominate the spectral energy distribution at low energies (~ 1 MeV). At 
high energies, the production spectrum can be significantly distorted by the effect of electromagnetic cascades. These effects 
are phase- dependent, and some variability modulated by the orbital period is predicted. 
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1. Introduction 

Many massive binary systems contain the final product of 
massive steUar evolution, a black hole or a neutron star. 
The mass transfer from the less evolved star onto the com- 
pact companion forms an accretion disk and generates X- 
ray emission. Some of these X-ray binaries also display 
coUimated outfiows of relativistic particles (jets) resem- 
bling scaled-down versions of the extragalactic quasars. It 
is widely accepted that these jets are powered by accre- 
tion processes and/or the rotational energy of the compact 
object. Their coUimation is thought to be related to the 
magnetic field structure in the surroundings of the com- 
pact object. 

Recently, microquasars (MQs) have been detected by 
ground-based Cherenkov telescopes at energies higher 
than 200 GeV (Aharonian et al. 2005; Albert et al. 2006: 
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Albert et al. 2007) . These detections have confirmed MQs 
as sites of effective acceleration of particles up to ultra- 
relativistic energies. The jets are the best candidates for 
generating the very high-energy radiation since they also 
present extended non thermal radio emission (see, e.g., 
Ribo 2005) and/or X-ray emission (see, e.g. Corbel et al. 
2002), evidencing the presence of a relativistic population 
of particles. Orbital modulation observed in some MQs at 
very different wavelengths and morphological changes re- 
veal a complex phenomenology, still not understood, tak- 
ing place in these gamma-ray binaries (e.g. LS I -1-61 303; 
Dhawan et al. 200fu) . For an updated review of the theo- 
retical aspects concerning MQs, see Bosch- Ramon (2007). 

The nature of the matter outflowing in the jets of MQs 
is still an issue of debate, although the presence of ions has 
been inferred from the detection of X-ray iron lines in the 



^ We note that these authors have proposed LS I-I-61 303 
as a pulsar binary system because the orbital variation of the 
extended radio emission morphology. For a critical assessment, 
see Romero et al. (2007). 
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case of SS 433 (Kotani et al. 1994; Migliari et al. 2002); 
hence, MQs jet models attempting to explain the broad- 
band spectral features of these sources have been centered 
in two different approaches. For a treatment considering 
the leptonic nature of the radiative processes that produce 
the primary emission, see e.g. Bosch- Ramon et al. (2006a). 
On the other hand, the so-called hadronic models, until 
now proposed only for high-mass systems, have predicted 
high-energy emission from the interaction between pro- 
tons accelerated in the jet and cold protons from the stel- 
lar wind of the primary star. The inelastic proton-proton 
collisions generate neutral pions that decay thereby pro- 
ducing the gamma-ray photons (see, e.g., Romero et al. 
2003; Orellana & Romero 2007). 

In hadronic models, the 7-ray emission should be ac- 
companied by a flux of high-energy neutrinos emerging 
from the decay of tt^. Their detection should be the 
unambiguous confirmation of the hadronic interactions 
(Aharonian et al. 2006). Also, since hadrons suffer energy 
losses to a lesser extent than leptons, the high-energy cut- 
off of the spectrum in hadronic models can be expected 
to occur at higher energies than in pure leptonic models 
(Orellana & Romero 2007). If there are other features that 
can be tested at lower energies to support the hadronic 
model over the pure leptonic one, they should be looked 
for in the secondary particle injection expected in such a 
scenario. 

We explore in this work the emission of leptonic ori- 
gin that is unavoidably produced in a hadronic scenario. 
There are two populations of secondary leptons that 
emerge in this context. On the one hand, relativistic lep- 
tons are injected along the jet by the charged pion decay 
channel. They are expected to be confined in the jet by 
the magnetic field and exposed to adiabatic and radia- 
tive energy losses of the same kind as in primary leptonic 
models. Regarding the radiative losses, we have accounted 
here for synchrotron energy losses under the magnetic field 
present in the jet and inverse Compton (IC) interactions 
with photons coming from the radiation field of the com- 
panion star. On the other hand, the substantial absorption 
of high-energy gamma-rays, mainly of neutral pion decay 
origin, in the photon field of the donor star injects en- 
ergetic pairs in the photosphere of the system. These 
leptons scatter stellar photons to high energies through IC 
interactions. The gamma-rays can be absorbed again, and 
an electromagnetic cascade develops in this way. 

The steady-state pair and photon distributions result- 
ing from the cascade reprocessing of the primary radiation 
can be estimated through numerical methods. We have 
performed Monte Carlo simulations of the electromagnetic 
cascade development along the line of sight. These one- 
dimensional cascades are induced by the primary hadronic 
emission when isotropization of very high-energy particles 
is not efficient and the magnetic field is rather weak, en- 
suring that the leptonic energy losses are IC-dominated. 
Variability of the outgoing spectra due to changes in the 
geometry of the IC scattering and photon-photon absorp- 
tion along the orbit, in which the inclination of the system 



plays an important role, can be determined. In what fol- 
lows we outline the model and present our calculations 
and results. 



2. The general picture 

The general picture treated in this work is the same as 
presented by Romero et al. (2003). The binary system 
is formed by a high-mass star and a compact object, 
here assumed to be a black hole. The star blows a radia- 
tively driven spherical wind that is partially accreted by 
the black hole and that developes an accretion disk. The 
mass accretion rate is assumed constant along the orbit, 
adopted here to be circular. Some of the accreted matter 
and energy is redirected into the outflowing jets. Hence, 
the total power of the jets is related to the total accretion 
energy budget. 

We assume a continuous jet, as observed in the sta- 
ble configuration of a MQ in a low-hard X-ray state (e.g. 
Fender et al. 2004). The launching mechanism for the jet is 
not treated here, and for simplicity, the jets are considered 
perpendicular to the orbital plane. These jets move away 
from the injection point at distance zq above and below 
the orbital plane and expand laterally as a cone with a jet 
semi-opening angle ~ 5°. They propagate with a mildly 
relativistic velocity characterized by a bulk Lorentz factor 
r = 1.5. Figure [T] shows a sketch of the configuration. 

The kinetic power of relativistic protons in each jet, 
is taken as a constant fraction (77) of the Eddington 
luminosity of the compact object. Two cases were con- 
sidered, labeled as Model A for 77 = 0.005 , i.e. L™^ = 
6 X 10^^ erg s^^, and Model B for 77 = 0.0005, i.e. 

= 6 X 10^^ erg s"^. We note that recently a low 
surface brightness arc of radio emission was discovered 
around the microquasar Cygnus X-1 and interpreted in 
terms of a shocked compressed shell of emitting gas in- 
flated by the jets (Gallo et al. 2005). This discovery led 
to the suggestion that the kinetic power of MQ jets can 
be as high as the bolometric X-ray luminosity of the sys- 
tem. For a system like Cygnus X-1, ^ 10^^ erg s~^; 
therefore, the values assumed above are consistent with 
this constraint. 

We assume that the wind ions penetrate the jet lat- 
erally, resulting in an effective production of charged and 
neutral pions through pp collisions. The main production 
region has roughly the size of the binary system, since 
farther away the wind dilutes and the number of inter- 
actions decreases rapidly. A penetration factor of 0.3 has 
been introduced to take the mixing between the jet and 
the wind material into account in a phenonienological way 
(see Romero et al. 2005). 

Concerning the wind that provides the proton targets, 
we consider a radiatively driven wind for simplicity. The 
radial dependence of the matter velocity is v{r) oc (1 — 
Rf/r) (Lamers & Cassinelli 1999). The particle density of 
this proton-dominated gas is obtained from the continuity 
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Table 1. Model parameters. The primes indicate quanti- 
ties measured in the jet co-moving frame. 





Jr^arameter, description lunitsi 


values 


e 


eccentricity 





a 


orbital semi-major axis [i?*] 


3 


i 


orbital inclination [°] 


30 


Mi, 


stellar mass [Mq] 


30 


M* 


stellar mass loss rate [Mq yr~"^] 


3 xlO"'' 






2 X 10* 


A/bh 


compact object mass [Mq] 


10 


Ro 


initial radius of the jet [-Rsch] 


5 


zq 


jet initial point [-Rsch] 


50 


X 


jet semi-opening angle tangent 


0.1 


a 


proton power-law index 


2 




proton minimum Lorentz factor 


2 




proton maximum Lorentz factor^ 


lO'^ 


r 


macroscopic Lorentz factor 


1.5 




stellar radius [Rq\ 


15 




stellar bolometric luminosity [erg s~^] 


6 xlO^* 




stellar surface temperature [K] 


3 X 10* 


U 


wind-jet penetration factor 


0.3 



2GA/bh/c2 



100 TeV 



equation. Along the jet (coincident with the z axis), it 
yields 



n(z) 



A'K'mp[z'^ + a?') 



1 - 



(1) 



The parameters adopted for the models are listed in Table 
[TJ Mostly, they correspond to a similar system to Cyg X- 
1: a high-mass binary with orbital period of 5.6 d, formed 
by an 09.7 lab super-giant primary (Walborn 1973) and a 
black hole of A^bh = 10 M©. The inclination of the orbital 
plane is taken to be i = 30°. The orbit is circular but its 
orientation induces periodic changes in the optical depth 
for photons produced close to the compact object. We 
follow the phase notation of Szotek & Zdziarski (2006): 
(/> = at the superior conjunction of the X-ray source, 
when the compact object is behind the primary. 

The jets are assumed to contain a randomly oriented 
magnetic field. We consider that, because of adiabatic 
expansion, the magnetic field follows i?jot ~ BQ{zf)/z), 
where zq = 50 i?sch- At the base of the jet we can con- 
strain the value of the magnetic field by its equipartition 
value, that is, B'^ ^ .^JsttU^, with the energy density 

C/p°' = Up^/{'!tRqc). Here, Ro is the initial radius of the 
jet, i.e. at zq- Then, B^'^ - 10^ G and B^"^ ~ 10^ G result 
for Models A and B, respectively. When emission from sec- 
ondaries in the jet is computed (see Section [5]), we adopt 
a magnetic field corresponding to 10% of the equiparti- 
tion values, incorporating an efhciency factor, as in mod- 
els with cold-matter dominated jets (Bosch- Ramon et al., 
2006a). 



3. Hadronic primary emission 

We outline here the basics of the hadronic model for pri- 
mary 7-ray production. The reader is referred to Romero 
et al. (2003 and 2005) for comprehensive descriptions. 

The particle spectrum of relativistic protons is as- 
sumed to be a power law in the jet reference frame with 
the canonical index (a = 2) from shock diffusive accelera- 
tion. The number density of particles entering the jet per 
unit of time can be determined as in Romero et al. (2003). 
The conservation of the number of particles makes the pro- 
ton flux depend on the distance to the compact object as 
Jp{z) oc {zq/z}'^. The Lorentz transformation that relates 
the co-moving jet frame with the observer's frame intro- 
duces an angular dependence in the proton injected flux 
(Purmohammad & Samimi, 2001). Speciflcally, we have 



Jp(Ep,d) = —[-l 



cKq/Zq^^ 

AnKz 



r 



-''+'(Ep^(3b^E^~mlc^ cost 
sin2 6l+r2fcos« 



where F is the jet Lorentz factor, and /3b the correspond- 
ing velocity in units of c. The constant Kq normalizes 
the energy distribution of the proton flux in the jet's co- 
moving frame. The zenith angle 9 is measured from the z 
axis. The gamma-ray photons are beamed along the pro- 
ton direction, reproducing their angular distribution. The 
maximum value of Ep obtained by balancing the accel- 
eration rate with the rate of energy loss by synchrotron, 
pp and p7 interactions exceeds the energy constraints im- 
posed by the equality between the maximum particle gy- 
roradius and the size of the accelerator, the latter taken 
as 0.1 times the jet radius at each height {Rjct{z) — x^)- 
In a model where the dependence of the magnetic fleld 
strength on z is inverse to that of Rjct{z), E^'^^ is con- 
stant along the jet, i.e. E^'''' = QAeB^Ro - 100 TeV. 
This value is used to impose an exponential cut-off on the 
above proton-flux. 

In our model, the generation of primary gamma-ray 
emission is not affected by the circular orbital motion 
of the compact object. With the parameters of Table [T] 
we have computed the spectral energy distribution (SED) 
of gamma-ray photons resulting from neutral pion de- 
cays. The (5-functional approximation is used to esti- 
mate the production rate of 7r°-mesonf@. Kelner et al. 
(2006) provide a recent and accurate approximation for 
the inelastic cross-section: ainc\{Ep) = (34.3 -1- 1.88 L-l- 



0.25 



1 - {EtijEp 



mb, 



where Et\, = 1.22 GeV 

is the threshold energy of production of Tr'^-mesons, and 
L = ln(£;p/lTeV). 

In Fig. [2] we show the hadronic gamma-ray luminosity 
generated along the jet (integrated along z) for different 
directions. The beaming factor is defined by the ratio 



^ For a comparison between the the 5-functional approxima- 
tion and other parametrizations of the 7-ray and emissivi- 
ties, see Domingo-Santamarfa & Torres (2005). 
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Fig. 1. Sketch of the configuration discussed in this paper. The main parameters are indicated. On the right a zoom-in 
of the region close to the black hole is shown. 

4. Secondary emission inside jets 

The expressions given by Kelner et al. (2006) are used 
to compute the leptonic energy distribution of injected 
along the jet by the charged pion decays. The resulting in- 
jection rate (leptons per units of energy, volume, and time) 
can be expressed in terms of the variable x = E^± / Ep as 




28 



60 120 180 

ep] 



9 10 11 

Log E [eV] 



12 



13 



14 



Fig. 2. Gamma-ray luminosity per unit of solid angle pro- 
duced by neutral pion decays in Model A (L™' = 6 x 10'^^ 
erg s^^). In decreasing order, the curves correspond to 9— 
0, 30, 60, and 90 degrees. The beaming factor is shown in 
the small panel. 



L{ff)/L{9 = 0) for energies greater than 10 GeV, where 
the angular dependence of the emerging flux is almost 
energy independent. 



<Ze(£^e) 



I'M 



dx 

,i{E,/x) Jp{Ejx) Fe{x, E,/x) —. (2) 







Here, Ff.{x, E^/x) is the spectrum of electrons from the 
t: fiiy^ decay, and it is described by 

C3(£;p)(lnx)2^3 



Fe{x, Ep) 
where 
Ci{Ep) = 

C2{Ep) = 
CsiEp) = 



Ci{Ep^ 



(1 

a;-h0.3a;(i-<^2(Bp)) 



1 



69.5 + 2.65 L 



0.3 L2 ' 
1 



(0.201 
0.279 -f 



- 0.062 L 
0.141 iH 



f 0.00042^2)1/4 
0.0172 2.2 



{~Hx)f, (3) 

(4) 

(5) 



0.3 -f (2.3 -f L)2 



(6) 



and L = ln(£;p/l TeV). 

For the injected proton flux (see Section 3), we obtain 
qe{Ee) 3 X 1027(i;e/eV)-i-92 leptons erg-^ a'^ cmT^, 
at the base of the jet, in the observer's reference frame 
for Model B (an order of magnitude higher for Model A). 
Once produced, we assume that the linear momentum of 
the leptons loses its angular dependence, erased by the 
random nature of the magnetic field. In the calculations, 
we integrate the 9 dependence of the differential rate of 
production of and treat the population as isotropic 
in the jet reference frame. The leptons injected along the 
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jet suffer from radiative cooling due to synchrotron and 
IC processes, as well as adiabatic expansion losses. The 
evolution of the particle energy distribution is studied by 
dividing the jet into slices of suitable size in order to con- 
sider the physical conditions as homogeneous within each 
of them. Afterward, to obtain the distribution evolution 
with we solve the continuity differential equation. Note 
that different evolved populations sum up at each height, 
since injection occurs all along the jet. The expressions 
used are given in Bordas et al. (2007). 

The secondary synchrotron emission is calculated in 
the presence of the magnetic field, which is determined by 
its value at the base of the jet B^. We compute the emis- 
sion for Bo = IQS G (Model A) and Bq = 10^ G (Model 
B). We assume the same lepton maximum energy, since for 
simplicity E'^'^^ is fixed to 100 TeV for both values of the 
magnetic field. For the specific emission and absorption 
coefficients, we used the expressions given by Pacholczyk 
(1970). We then computed the spectral energy distribu- 
tions and transformed them into the observer frame. 

As a first-order approach, we assumed that the IC 
takes place for any lepton energy in the Thomson 
regime and under an isotropic photon field. Actually, the 
Thomson approximation is not valid for particles with en- 
ergies above ^ 10-100 GeV, i.e. at the highest energies of 
the secondary leptons produced in this scenario. Above 
these energies, the IC interactions occur in the Klein 
Nishina (KN) regime. This fact does not affect the radi- 
ated synchrotron spectrum significantly, since it is mainly 
produced in the inner regions of the jet, close to the com- 
pact object, where the synchrotron process is the main 
cooling mechanism; i.e. the particle spectrum is not notice- 
ably affected by KN cooling there. The KN effect would 
otherwise affect the highest energies of the IC spectrum, 
but this energy range of the spectrum is dominated by 
neutral pion decay 7-ray photons (see next section). 

In Figure [3l we show the synchrotron and IC spec- 
tra. These calculations allow us to estimate the power go- 
ing to this radiative channel in order to compare it with 
other components like IC emission from electromagnetic 
cascades or synchrotron emission from the secondaries in- 
side the jet. The set of explored parameters illustrate the 
impact of the IC/synchrotron energy loss balance in the 
secondary broadband spectrum. It is seen that the syn- 
chrotron luminosity is at most ~ lO'^^ erg s~^, and well 
below the IC emission for the both values of the adopted 
magnetic field. 

5. IC cascades along the line of sight 

The radiation field of early-type stars provides a suitable 
target for the absorption of gamma-ray photons. In Fig. H] 
we show the optical depth as a function of the gamma- 
ray energy for the injection at height z — a/2. Hereafter 
we consider that as a representative point for the injected 
primary emission. The plotted optical depths were calcu- 
lated as in Dubus (2006). The absorption probability is 
modulated on the orbital period of the binary as a result 
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- . Model A, B(,=10'^G X 
■ Model B,Bo=10''g / \ 




Log E [eV] 

Fig. 3. Synchrotron and inverse Compton emission from 
the leptons produced inside the jet by the charged pious' 
decay. The inverse Compton interactions are calculated in 
the Thomson regime. 

of the phase dependence of the geometric path (see the 
sketch in Szostek & Zdziarski 2006). 

Energetic electron-positron pairs are materialized by 
photon-photon interactions. These leptons, in turn, boost 
the stellar photons to high-energies via IC scattering. The 
7-ray absorption and production mechanisms can proceed 
very fast, resulting in the development of an electromag- 
netic cascade (Bcdnarek 1997, 2000, 2006; Aharonian et al. 
2006). As a result, the energy of the original photons is 
distributed between a certain number of secondary parti- 
cles and photons with lower energy. 

We computed the 7-ray spectra formed in cascades 
traversing the anisotropic stellar radiation field. Other ra- 
diative fields (like the accretion disk field) are less im- 
portant in the present context. We implemented a Monte 
Carlo simulation code based on the scheme outlined by 
Protheroe (1986) and Protheroe et al. (1992). A descrip- 
tion of the treatment is given in Orellana et al. (2006). 

The injected flux that initiates the cascades is obtained 
by fitting the gamma-ray spectrum for i = 9 = SO de- 
grees. For Model A, we fLiiSlog L{E-f) [ergs~i] - 32.27-1- 
0.06 log £;t, [eV] for 1 GeV < < 10 TeV. The absorb- 
ing radiation field has a black-body spectrum. It comes 
from the 0-type star, characterized by Tcff ~ 30000 K 
and ~ lO^^ ^m. 

For simplicity, any effects of the presence of a mag- 
netic field are disregarded in these simulations. However, 
we note that the cooling rates by synchrotron and IC pro- 
cesses are equal at a small value of the magnetic field, 
namely i? ~ 0.1 — 1 G for electron energies ~ 1 TeV, deep 
in the KN regime. Hence, the value assumed for the mag- 
netic field is crucial for the development of the cascade. 

The isotropic luminosity is quoted here, being defined 
L(E) =A-KdL{E)/dn. 
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Fig. 4. Optical depth for gamma-ray photons injected at 
z = a/2. The absorber photon field generated by the donor 
star is anisotropic, considered here as a black body emit- 
ter. Each curve corresponds to a fixed orbital phase, with 
= at the superior conjunction of the X-ray source, 
when the compact object is behind the primary. 



The geometry and radial dependence of the stellar mag- 
netic field close to the star has been derived by several 
authors (e.g., Eichler & Usov 1993) under different sim- 
plifying assumptions. The main parameter for B{r) is the 
value of magnetic field at the stellar surface, which is not 
well known, but it is thought that it could reach ~ 10^ G 
at most. Then, the field could be in the range 10~^ — 10 G 
in a region at a few stellar radii around a typical O stai0. 

As the cascade evolves, departing from the injection 
point toward the observer, it will certainly reach regions 
where the stellar magnetic field has decayed; hence, there 
is still some room for our assumption (small or negligible 
B) to be valid. Effects of the magnetic field influencing the 
pure IC cascade by changing the lepton direction of prop- 
agation have been incorporated by Bednarek (1997 and 
subsequent works) to cascade simulations. We neverthe- 
less assume that, at the lepton energies to which most of 
the cascade develops, particles radiate before isotropizing 
after their creation. 

The farthest distance from the source up to which we 
follow the cascades is ~ 10 times the orbital separation. 
Photons that pass beyond that distance are stored to form 
the outgoing spectra. At distances greater than the or- 
bital separation, the cascade efficiency is reduced by the 
dilution of the stellar radiation field. The spectra are ob- 
tained by constructing histograms f{E), sorting photons 
into bins of width A(log£') = 0.5. Some emerging spectra 
are shown in Fig. [5l 

The injected primary power is clearly re-distributed 
to lower energies. At energies higher than ~ 10 GeV, the 
spectra do not follow a simple power-law, but present a 



* These values of B axe taken from Benaglia & Romero 
(2003) for colliding wind regions of WR-I-OB systems. 
Concerning the stellar magnetic field the physical context is 
the same than the one treated here. 
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Fig. 5. Isotropic luminosity of cascade spectra for relevant 
orbital phases. The injected (hadronic) primary gamma- 
ray power is indicated. Thick lines correspond to Model 
A and the thin ones to Model B. For the orbital phases, 
= when the compact object is behind the star, (j) = 0.5 
when it is in the front. 



more complex structure, much softer than the injected 
spectrum. At energies below ^ 100 MeV, the cascade ef- 
fect is essential for raising the luminosity, and the phase 
dependence is stronger than at higher energies. The results 
point to a time dependence anti-correlated for ~GeV sig- 
nal compared with that at ~ TeV. 



6. Discussion and summary 

We have considered the simplified geometry of a circu- 
lar orbit and a spherical wind for a high-mass MQ in 
order to illustrate the secondary emission produced in a 
hadronic model. We considered a power-law distribution 
with a spectral index a = 2 for the relativistic proton 
flux injected at the base of the jet. The gamma-ray and 
leptonic emissivities from 7r-meson decays were computed 
using the latest available analytic approximations for the 
inelastic pp cross-section (Kelner et al. 2006). Relavistic 
leptons that are produced inside the jet via charged pion 
decays radiate subsequently via synchrotron and IC pro- 
cesses. 

The primary hadronic emission is found to be strongly 
affected by the pure IC cascade development in the photo- 
sphere of the primary star. We note that the synchrotron 
dominance of the leptonic energy losses can prevent the 
pure IC cascade from developing. In our treatment, we re- 
strict ourselves to a low magnetic field in the region where 
the cascade occurs. In the case of a strong magnetic field, 
cooling via synchrotron radiation of the leptons injected 
by 7 — 7 absorption will suppress the cascade development 
(e.g. Khangulyan et al. 2007). 

The treatment implemented for the simulation of the 
electromagnetic cascades ignores the effect of any other 
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process than IC interactions in the direction to the ob- 
server. This is a good approximation for when particles 
radiate before isotropize. Efficient isotropization of ~ 0.1- 
10 TeV particles would require diffusion coefficients close 
to the Bohm one and magnetic field values 1 G. In 
addition, Doppler boosting can also favor the jet radia- 
tion to dominate. Nevertheless, low-energy particles ab- 
sorbed in different directions from that of the observer 
can present lower energy emission re-directed toward the 
observer since their direction of motion could indeed be 
randomized. Approximate estimates of the emitted radi- 
ation by leptons confined relatively close to the binary 
system are given in Bosch-Ramon (2007). 

The cascade reprocessing can lead to significant emis- 
sion at energies well below those originally injected (with 
primary energy > 1 GeV, due to the pion-creation 
threshold). The total radiation produced inside and out- 
side the jet by charged pion and cascade secondaries is 
similar to that produced at gamma rays by the neutral 
pion decay. Dctectability under the parameters assumed 
here can be discussed at different bands. 

At radio frequencies, the presented spectra are signif- 
icantly hard, yielding a flux ^mJy at distances of a few 
kpc, which represents ~ 10% of the fluxes found in MQs 
in the low-hard state, i.e. 15 mJy (for Cygnus X-1, 
see, e.g., Pooley ct al. 1999). If this radio emission pro- 
duced by secondaries could be distinguished when ob- 
served at very high resolution, it would appear morpho- 
logically as a jet-likc source surroimdcd by a radio halo, 
because of the highly beamed relativistic leptons gener- 
ated by charged pion decays, on the one hand, and those 
generated by pair creation in the stellar photon field on 
the other. From these points, it seems that a primary lep- 
tonic component would be needed to explain the observed 
radio fluxes. Nevertheless, electrons have lifetimes that are 
long enough to reach regions of more suitable conditions 
to emit in the radio band, and the secondaries created 
in the jet could also be rcaccelerated once they have lost 
their energy. Thus, it cannot be ruled out that extended 
radio emission is in fact produced by secondaries, in which 
case a primary electron population might not be necessary 
for explaining the non thermal spectrum observed in these 
systems. Such a case is outside of the assumptions made 
along this work, but could be interesting to explore in the 
future. 

At soft X-rays, synchrotron emission from secondaries 
could be detectable by the current X-ray instruments, al- 
though low-hard state MQs used to present significant 
X-ray emission produced in the accretion disk and the 
corona, which could mask this non thermal contribution 
to some extent. 

At hard X-rays, it is still not clear whether the radia- 
tion observed in some objects has a thermal or non ther- 
mal origin (see, e.g. Markoff et al. 2005; Maccarone 2005; 
and also, Bosch-Ramon et al. 2006c). In any case, the lu- 
minosities at hard X-rays obtained here by summing dif- 
ferent contributions would not suffice for explaining those 
observed in some MQs during the low-hard state. 



The radiation in the MeV range seems to be non neg- 
ligible, reaching ^ 10'''^ erg s"-'^, although it still appears 
too low to be detected on a reasonable time span for the 
present instruments. At high and very high-energy gamma 
rays, GLAST and modern Chcrcnkov telescopes like HESS 
or MAGIC could detect the source, which could present a 
relatively soft spectrum due to electromagnetic cascading. 
In fact, while finalizing the research for this paper, the de- 
tection of Cygnus X-1 by the MAGIC collaboration was 
announced (Albert et al. 2007). Throughout this work, it 
was not our purpose to explain the detailed behavior of 
this specific source; instead, we used some of the param- 
eters of Cyg X-1 to apply the general treatment to this 
very interesting example. However, it is remarkable that 
the source has been detected by MAGIC around phase 0.9, 
when photon-photon absorption is high. Actually, it may 
be related to some additional variability of the production 
mechanism not accounted for here but linked to this phase 
(e.g. aff'ected by perturbations in the stellar wind mate- 
rial; e.g. Feng & Cui 2002), although the investigation of 
such a possibility is left for future work. 

Even from the simple model considered here, we have 
been able to discuss variable emission in some energy 
bands. Many specific factors should lead to more complex 
variations (see e.g. Bosch- Ramon et al. 2006a). For in- 
stance, realistic modeling of the variability would include 
a variable accretion rate (e.g. induced through the vari- 
able nature of the outflows of the donor star, or simply 
through an eccentric orbit) and possible changes through 
the misalignment between the jet and the orbital axis (i.e. 
a precessing jet, see Kaufman Bernado et al. 2002). In 
such a more complex scenario, it woiild also be worthwhile 
study the evolution along the orbit of the IC hard X-rays 
and soft gamma rays due to the angular dependence in 
the cross section of this mechanism (see Khangulyan et 
al. 2007 and references therein). 
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